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Abstract 

This Letter investigates the compressible energy transfer process. We extend a 
methodology developed originally for incompressible turbulence and use databases 
from numerical simulations of a weak compressible turbulence based on Eddy-Damped- 
Quasi-Normal-Markovian (EDQNM) closure. In order to analyze the compressible 
mode directly, the well known Helmholtz decomposition is used. While the compress- 
ible component has very little influence on the solenoidal part, we found that almost 
all of compressible turbulence energy is received from its solenoidal counterpart. We 
focus on the most fundamental building block of the energy transfer process, the 
triadic interactions. This analysis leads us to conclude that, at low turbulent Mach 
number, the compressible energy transfer process is dominated by a local radiative 
transfer (absorption) in both inertial and energy containing ranges. 
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1 Introduction 


Compressible turbulence research, especially using direct numerical simulations 1-5 , 
was invigorated recently by its engineering applications such as the high speed civil 
transport and supersonic combustions ram-jet engines. It is well known that direct nu- 
merical simulation is restricted to very low Reynolds numbers. The Reynolds average 
and large-eddy simulation are two alternative methods important for high Reynolds 
number scientific and engineering calculations. The accuracy of these computations 
relies on their respective turbulence models. Model development, in turn, requires a 
fundamental understanding of the energy transfer process. This Letter reports a first 
in-depth investigation of the energy transfer process in homogeneous isotropic com- 
pressible turbulence. We expect that the results of this work will find use in future 
model developments for Reynolds averaged models and large-eddy simulations. 

2 Analysis 

To analyze the compressibility effects, we use the Helmholtz decomposition to split 
the velocity vector into a solenoidal part u 5 (K,t), which corresponds to the velocity 
fluctuations perpendicular to the wave vector K in the Fourier space, and a com- 
pressible part u c ( K, t), which corresponds to fluctuations in the direction of the wave 
vector. We assume that the fluid is barotropic and that compressibility of the turbu- 
lence is weak. Moyal 6 pointed out that in the absence of mean velocity gradient, the 
interactions between these two components are exclusively due to the nonlinear terms 
and are very important. To obtain high Reynolds number flow fields, we shall use the 
energy transfer equations for weak compressible turbulence based on the EDQNM 
closure 7-9 , a two-point statistical theory. A forcing is applied in the large scales of 
the solenoidal velocity component. We extend the methodology developed originally 
for incompressible turbulence 10-16 and carry out our analysis when both solenoidal 
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and compressible modes have reached their asymptotic states 17 . 

In the framework of our hypothesis, the governing evolution equations for a 
weakly compressible turbulence (M t <C 1) are as follows: 

(1) The spectrum ( E ss ) of the “solenoidal” velocity correlation: 

E SS (I< } t) = -2 vK 2 E ss (K, t) + T SS (I< } t). (1) 

(2) the spectrum ( E cc ) of the “compressible” velocity correlation: 

^ E CC (K , t) = -2EK 2 E cc (K } t) + T CC (K , t) - E CP (K } t). (2) 

Here < p > is the mean density. In the case of a Stokes fluid, v' = yrjyy = | zy 
where p and A are the two dynamic viscosities. The evolution equation for E cc 
requires the velocity-pressure correlation term, E cp (K,t ), which is given by a sepa- 
rate transport equation that depends on the potential energy. This term, along with 
the viscous terms, transfers the compressible kinetic energy into internal energy. We 
will not investigate it here since it is not relevant to the nonlinear energy transfer of 
compressible turbulence. 

In this letter, we report only results for a turbulent Mach number of 10 -2 and 
a Taylor micro-scale Reynolds number of 140. The turbulent Mach number, M t = 
VF/ c o, is a dimensionless number characterizing the effects of compressibility and 
the Reynolds number is defined as R e = where q 2 is twice the turbulent kinetic 

energy per unit of mass, c 0 the sound speed and e the dissipation. In the next figures, 
the spectra of double correlations will be given in m 3 s~ 2 and K is in to -1 . 

Figure 1 demonstrates the existence of an extended solenoidal inertial range 
spectrum, the same as that of Kolmogorov for incompressible turbulence. In the 
wave-number range corresponding to the inertial range of the solenoidal mode, the 
compressible spectrum shows a slope 17 of —11/3. 

In the equations (1) and (2), T ss (k ) and T cc (k) are the nonlinear energy transfer 
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terms and they are composed of different contributions (details given in Appendix): 

r ss (*0 = XA ss «. (3) 

i=l 

T cc (k) = J2Tf c (k). (4) 

8 = 1 

It is interesting to note that Tg S (k ) = Tf s (k) + T^ s (k) is the contribution from the 
purely solenoidal velocity component and is the same as the incompressible turbu- 
lence transfer function. In fact, T ss is dominated by Tg S for turbulent flows at low 
Mach numbers and it is the same as that found in incompressible studies 10-16 (not 
shown). Other terms in T ss (k ) are the energy transfer resulting from the interac- 
tions between the solenoidal and compressible modes. Therefore, the remaining terms 
in T ss , denoted as Tq S hereafter, are negligible when one examines the solenoidal 
velocity transfer 8,18 . 

We now focus on the energy transfer process in T cc (k), which is a pure compress- 
ible transfer term. We stress that Tq S term plays important and fundamental roles 
in exchanging energy between the solenoidal and compressible modes when our fo- 
cus is the compressible energy transfer. Indeed, the compressible component energy, 
T CC (K), has the same magnitude, but opposite sign from that of Tq S (K). Note that 
the compressible transfer function is positive for all resolved wavenumbers (Figure 
2). This, in turn, indicates that all of the compressible energy is transfered from 
the solenoidal mode. To distinguish it from the well known “cascade” picture and in 
analogy with the radiation transfer, we shall call the positive peak absorption while 
the negative peak is called emission. 

The most fundamental building block of energy transfer process is the triadic 
interactions. Specifically, we are interested in the energy transfer for a given mode 
K due to its interactions with all the pairs of modes P and Q = K - P that form 
a triangle with K. For this reason, we introduce the triadic energy transfer function, 
T(K, P, Q), according to 


3 



( 5 ) 


T SS (K)= J2 T ss (K,P,Q). 

P,Q = |K-P| 

T CC (K)= J2 T cc (K,P,Q). (6) 

P,Q = |K-P| 

Here T(K } P } Q) is defined as energy transfer to K due to triads with one leg in 
Q and the other in P. An examination of the purely incompressible contributions 
(Tg S (K, P, Q) ) reproduces the results of incompressible turbulence 10-16 and indicates 
again that the solenoidal triadic energy transfer is not affected by compressible ef- 
fects. The total triadic solenoidal transfer T SS (K , P, Q) (with the compressible terms 
Tq S (K, P, Q) included) is essentially the same 19 . As a result, the compressibility has 
very little influence on the solenoidal triadic interactions. 

We now turn our attention to the triadic interactions in compressible energy trans- 
fer, T cc (K, P, Q). In Figure 3, we present T cc (K, P, Q ) for various Q values when P 
is in the inertial range (P = 512). We first observe that the structures of T cc (K, P, Q ) 
are rather similar for different Q values. All of them show the radiative (absorption) 
type of energy transfer. Figure 4 illustrates the same kind of behavior for P in the 
energy containing range (P = 128). Recall that there is only a direct energy cascade 
in incompressible energy transfer. 

We now examine the behavior of the function T cc (K, P) which is given by 

T CC (K, P)=J2 T CC (K, P, Q). 

Q 

Here T CC (K,P) is defined as the transfer of compressible energy to wavenumber K 
due to triads with at least one leg in a wavenumber P. Of particular importance, 
it indicates the direction as well as the locality of the energy transfer. First of all, 
T cc (K, P) is dominated by energy absorption when P is either in the energy con- 
taining or in the inertial range. Since most of the energy transfer is concentrated 
about P = K (within a decade), we conclude that the compressible energy transfer 
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is relatively local. We note that there is no further cancelation when summing over 
all Q in T CC (K, P, Q). Moreover, at this low turbulent Mach number, the compress- 
ible energy input is dominated by the radiative energy transfer from the solenoidal 
modes; there is no significant compressible energy cascade from large scales to the 
small scales. The dominance of energy absorption clearly establishes the different 
physical mechanism of compressible energy transfer compared to its incompressible 
counterpart. 

3 Conclusions 

This letter focuses on the influence of the compressibility effects on compressible 
energy transfer. Another important issues, the detailed analysis of the individual 
terms in Eqs. (5) and (6) as well as the locality of the interacting scales, are out of 
scope of present letter and the reader is referred to our full paper 19 for details. At 
low turbulent Mach number, our observation suggests that the compressible energy 
transfer process is completely different from its incompressible counterpart. The 
well known “cascade” picture of the incompressible local energy transfer is no longer 
dominant here. We believe that this result should be reflected in all models for 
compressible turbulence. 
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APPENDIX 


The transfer terms T ss and T cc are given by : 


rj-iSS rj-iSS _|_ rj-iSS _|_ rj~iSS _|_ rj~iSS _|_ rj~iSS 


Jlbb Jlbb | Jlbb | JlOO | JlOO | JlOO | JlOO 


( 7 ) 

( 8 ) 


The different contributions to T ss are : 

if* -Ijq 1 ~ XyZ ~ 2yV el^- SS E ss (PJ)E ss (,Q,t)dPdQ (9) 

Ti s =J^ (1 -^ + !/ 2) offf q ss - cc E ss IPJ)E cc I.QJ)dPdQ (10) 

T| s = -f^( xy + z 3 )4lf h ? s - ss E ss (K.t)E ss (Q.t)dPdQ (11) 

Tf = / a ^ {2xy)»f,^ c - ss E ss IJ<J)E cc l,QJ)dPdQ (12) 

If* ~~!j^ W1 - z 2 ))0pqJ/ s ~ ss E ss (K,t)E ss (Q,t)dPdQ (1:1) 
The different contributions of T cc are : 

7f c (A',i) = ^ tff ((* + »*) , )»g55**-**£!**(P,«)£!**(<?,«)<tf>d<? (14) 

7? C (A',() = i % TT^vf ^ S - CC E SS (P.t)E cc (Q,t)dPdQ (15) 
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(16) 


T. 


cc 


r K 3 

<*•*> = Lpq {: 


2 \nCC-CC-CC 
rKPQ 


E cc (P,t)E cc (Q,t)dPdQ 


Tf°(K,t) = -/ a P 2*(1 - .- 2 )^R SS - cc £ cc (A',i)£ ss (Q.i)<lPflQ (17) 


lf c (K,t) — (2- 3 - - + xy)6%]f s - cc E ss (Q,t)E cc (K,t)dPdQ (18) 

Tf c (K,t) = J a E (2xy)0^ c - cc E cc (K,t)E cc (Q,t)dPdQ (19) 

The integration in the P, Q plane extends over a domain such that K, P and Q 
can be the three sides of a triangle. 

In the equations, the polynomial expressions in x, y, z are coefficients associated 
with the geometry of the triad, x, y, and z are the cosines of the angles respectively 
opposite to K, P, Q in the triad K, P and Q. 9 are the decorrelation times involved 
in the model. Their expressions can be found in Bataille [8]. 

We should remark that the transfer equations are not completely given since we 
have omitted the E cp contributions. The complete equations can be found in Bataille 
[ 8 ], 
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